Chaos Theory 


Unpredictability of nonlinear systems 


, Chaos ( : Χάος, "Chasm")!"! was the first thing to exist: "at first Chaos came to be" (or was)!” "but next" (possibly out of Chaos) came 
,and (elsewhere the son of ). Unambiguously born "from Chaos" were (Darkness) and (Night). 


(a) Response of the annual mean surface temperature and (b) precipitation to 
Special Report on Emission Scenarios A1B emissions, in 21 climate models 
that contributed to the IPCC Fourth Assessment Report. 
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Schematic of a probabilistic weather forecast using initial condition 
uncertainties. 
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Example of 66h probabilistic forecast for 15-16 October 1987. 
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Two contrasting ensemble seasonal forecasts from the European Centre for 
Medium-range Weather Forecasts (ECMWF) for the evolution of El Nino. 


a b 
(a) NINO3.4 SST anomaly plume 0) NINO3.4 SST anomaly plume 
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Laplace's demon 


The demon of Laplace (1794-1827) 


We may regard the present state of the universe 
as the effect of its past and the cause of its 
future. An intellect which at a certain moment 
would know all forces that set nature in motion, 
and all positions of all items of which nature is 
composed, if this intellect were also vast enough 
to submit these data to analysis, it would 


embrace in a single formula the movements of 
the greatest bodies of the universe and those of 
the tiniest atom; for such an intellect nothing 
would be uncertain and the future just like the 
past would be present before its eyes. 

A Philosophical Essay on Probabilities, 1814 
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APPLICABILITY OF 


these four initial values, the trajectory can be integrated 
to the next point satisfying (9), which is P.] It can 
also be shown that the mapping is area-preserving [ see, 
e.g., Birkhoff (1927, p. 152); and see Moser (1962) for 
an important theorem concerning such mappings |. 


3. RESULTS 


After some trials, the following potential was chosen 
for study: 


)= 3 (a4- yt-2x2y — 35) (11) 


because: (1) it is analytically simple; this makes the 
computation of the trajectory easy; (2) at the same 
time, it is sufficiently complicated to give trajectories 
which are far from trivial, as will be seen below. It seems 
probable that the potential (11) is a typical representa- 
tive of the general case, and that nothing would be 
fundamentally changed by the addition of higher-order 
terms. 


Y| c+ 0.0833 


Γιο. 4. Results for E = 0.08333. 


Figure 2 shows the equipotential lines, Near the 
center they tend to be circles; farther out they become 
elongated in three directions. The particular equi- 
potential U=} consists of three straight lines, forming 
an equilateral triangle. 

A number of orbits were computed by numerical 
integration of the equations of motion: 


і---90/Әх---х--2ху, 
й--9П/ду--у-х--у, (12) 


Ав а check, some of the orbits were computed inde- 
pendently by each of us, using different computers 
(CDC 1604 and IBM 7090) and different integration 
schemes (Adams and Runge-Kutta). The following 


results were obtained using the Runge-Kutta method ; 
during the numerical intecratinn the enerov wac nh- 
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Fic. 5. Results for E=0.12500. 


Figure 3 shows a set of points P; for a typical tra- 
jectory. They seem to lie exactly on a curve. In fact, 
more points have been computed than those plotted 
here; after the 150th point there is still no perceptible 
deviation from a curve. It may be interesting to remark 
that the successive points Py, Р», Pi... (represented 
here by 1, 2, 3---) rotate regularly around the curve. 
Тһе figure is topologically identical to one where the 
points P; would lie on a circle of center O, the angle 
between OP; and ОР; having a constant value a. This 
constant is not the same for different trajectories. In the 
case of Fig. 3, its approximate value isa = 0.1143 (taking 
one revolution as the unit). « is generally not rational, 
so that no point P; will come back exactly on the initial 
point Рі, and the infinite set of the points P; is dense 
everywhere on the curve. Па happens to be a rational 
number ῥ/ῃ, the point Pei will be identical with Р; and 
the orbit is periodic. 

Figure 4 shows the complete picture in the (у,)) 
plane, for a given value of the energy: / = yy = 0.08333. 
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Modeling of the weather 


& Edward Lorenz (1962) 
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Butterfly Effect in the Lorenz 
System 


Examples of finite-time error growth оп the Lorenz attractor for three 
probabilistic predictions starting from different points on the attractor. 
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Figure 12. Cobweb plots of the logistic map pulling initial population values of 0.1 (A), 0.5 (В) and 
0.9 (C) into the same fixed-point attractor over time. At this growth rate parameter value of 2, 
the Lyapunov exponent is negative. 
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Figure 6. Phase diagrams of the logistic map over 200 generations for growth rate parameter values of 
2.9 (A); 3.5 (B); 3.56 (C); and 3.57 (D). When the parameter is set to 2.9, the model converges at a single 
fixed-point. When the parameter is set to 3.5 or higher, the model oscillates over four points, then eight, 
and on and on as it bifurcates. 
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Figure 10. Phase diagrams of the time series in Figure 9. (A) is a two-dimensional state space version 
(the xy-plane) of the three-dimensional one (B). 
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Figure 9. Plot of two time series, one deterministic/ chaotic from the logistic map (blue), and one 
random (red). 
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Figure 1. Time series graph of the logistic map with seven growth rate parameter values over 


20 generations. 
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Figure 2. Bifurcation diagram of 100 generations of the logistic map for 1000 growth rate parameter 
values between zero and four. The vertical slice above each growth rate depicts the system's attractor 


at that rate. 
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Figure 3. Bifurcation diagram of 100 generations of the logistic map for 1000 growth rate parameter 
values between 2.8 and 4. The vertical slice above each growth rate depicts the system's attractor at 
that rate. 
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Figure 4. Bifurcation diagram of 100 generations of the logistic map for 1000 growth rate parameter 
values between 3.7 and 3.9. The system moves from order to chaos and back again as the growth rate 


is adjusted. 
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Figure 5. Bifurcation diagram of 100 generations of the logistic map for 1000 growth rate parameter 
values between 3.84 and 3.856. This is the same structure that we saw earlier at the macro-level in 


Figure 3, because chaotic systems' strange attractors are fractal. 
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Figure 7. Cropped phase diagrams of the logistic map over 200 generations for: (A) a growth rate 


parameter value of 3.9; and (B) 50 growth rate parameter values between 3.6 and 4 (the chaotic regime), 
each with its own colored line 
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Figure 8. Cobweb plots of the logistic map for growth rate parameter values of: (A) 1; (B) 2.7; (C) 3.5; 
(D) 3.9. The diagonal gray identity line represents y = x; the red curve represents the logistic map 
as y — f(x) for each of the four parameter values; and the blue cobweb line represents the system's 
trajectory over 100 generations. 
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Figure 13. Plot of two time series with identical dynamics, one starting at an initial population 
value of 0.5 (blue) and the other starting at 0.50001 (red). At this growth rate parameter value of 
3.9, the Lyapunov exponent is positive; thus, the system is chaotic, and we can see the nearby points 


diverge over time. 
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The Lorenz Attractor 
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Figure 9.5.2 


At first the trajectory seems to be tracing out a strange attractor, but eventually it stays 
on the right and spirals down toward the stable fixed point С”, (Recall that both С” 
and С” are still stable at r = 21.) The time series of y vs. t shows the same result: an 
initially erratic solution ultimately damps down to equilibrium (Figure 9.5.3). 


Figure 9.5.3 
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This solution predicts that the waterwheel should ultimately rock back and forth 
like a pendulum, turning once to the right, then back to the left, and so on. This is 
observed experimentally. m 
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Attractor reconstruction 


Roux et al. (1983) exploited a surprising data-analysis technique, now known as 
attractor reconstruction (Packard et al. 1980, Takens 1981). The claim is that for 
systems governed by an attractor, the dynamics in the full phase space can be re- 
constructed from measurements of just a single time series! Somehow that single 
variable carries sufficient information about all the others. 

The method is based on time delays. For instance, define а two-dimensional 
vector x(t) = (Bit), Bit + 1)) for some delay т > 0. Then the time series B(r) gen- 
erates a trajectory x(r) in a two-dimensional phase space. Figure 12.4.2 shows the 
result of this procedure when applied to the data of Figure 12.4.1, using T = 8.8 
seconds. The experimental data trace out a strange attractor that looks remarkably 

like the Rössler attractor! 
Roux et al. (1983) also considered 
the attractor in three dimensions, by 


defining the three-dimensional vec- 
tor x(t) = (Bit), Ва + 7), Ва  21)) 
To obtain a Poincaré section of the 


attractor, they computed the intersec- 
tions of the orbits x(r) with a fixed 
plane approximately normal to the 
orbits (shown in projection as a 
dashed line in Figure 12.4.2). Within 
the experimental resolution, the data 
fall on a one-dimensional curve. 
Bit) Hence the chaotic trajectories are 
confined to an approximately two- 


Figure 12.4.2 Roux et οἱ. (1983), p. 262 dimensional sheet. 
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1 The S&P 500 long term chart on 35 years shows the triple top 300 
breakout in 2013. That was an incredibly important event in the 
stock market history. Based on this view, the S&P 500 has much 
more upside potential. However, it is mandatory to always check a 
** chart on a larger timeframe which can reveal the bigger picture, 
i.e. more important trends, trend lines, patterns. 
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Figure 6 


Metabolic rate for a series of mammals and birds as a function of mass. The scale is logarithmic, and 
the slope of 3/4 exemplifies Kleiber's law: The metabolic rate of an animal of mass m is proportional 
to m. This law has recently been explained by West et al. 1997. Figure taken from West et al. 2000. 
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Admin. services: B = 1.08(0.02), In(Y,) = -3.61(0.23) 
Wholesale brokers: B = 1.29(0.03), n = -8.91(0.39) 
= Creatives: B = 1.11(0.01), In(Y,) = -3.23(0.14) 
Inventors: $ = 1.47(0.06), TAA = -14.09(0.74) 
= Burglary: B = 1.01(0.02), in(Y,) = -5.16(0.24) 
Robbery: В = 1.35(0.03), In(Y,) = -11.65(0.34) 
«= High school: Ё = 100), ПА) = -0.15(0.04) 
Graduate: $ = 11100.02), А = -3.82(0.21) 
= Chlamydia: f = 1.06(0.02), п(Ү,) = -6.37(0.29) 
Syphilis: $ = 1.46(0.05), In(Y,) = -16.91(0.65) 


Scaling and Politics 


* Democrats, B=1,13720.004%, R2-0.98 
a Republicans, B=0.916+0.007, Е2-0.96 


z 


5 
> 
5 
+ 


Scaling in Complex 
etworks 


Degree distribution of WWW 


Degree distribution of Internet 
routers 


Network Network 
172.16.0.0 


Network Destination Exit 
Protocol Network Interface 
EIGRP 10.1.1.0 ҒА0/1 


OSPF 172.16.0.0 ΕΑΟ/2 


Degree distribution in Social 
Networks 


_ Slope: -1.514 


Slope: -1.611 


Slope: -2.602 


Gibrat’s law 


® Robert Gibrat (1931) 


® Growth rate (x=company- asset- or city size) 


© 
2 
— 
= 
Ξ 
5 
5 


Fiscal year 


Growth is independent of 
size 


Probability density 
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Log growth rate r Log growth rate r 
& Under this assumption, size distribution is power 
law again with a approximately 1. 
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Frequency 
© 
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10 ю 


| ! 1 
10: 10: 10% 
Firm size (employees) 


Log frequency Inf(S) versus log size In$ of U.S. firm sizes (by number of employees) for 1997. 
Ordinary-least-squares fit gives a slope of 2.059 (s.e. — 0.054; R? — 0.992). This corresponds to a 
frequency f (S) ~ $~*°*, i.e., a power law distribution with exponent ζ = 1.059. This is very close to 
Zipf’s law, which says that ( = 1. Figure taken from Axtell 2001. 
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(а) Sales 


n — 456,712 
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Firm size ( Y) 
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Тһе Black Swan 


Nassim Taleb 


Fat Tail 


Normal Distribution 
(Bell Curve) BEN C = Standard Deviation 


Fat Tail 
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Fig. 2. Earthquake magnitude distribution showing a power-law behavior 
over 6 decades. The graph follows logıo N(M > m) « — bm, where b is the 
Gutenberg-Richter exponent b = 1 (dashed red line). The roll-off for m < 2 is 
due to difficulties with detecting very small earthquakes. 
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$ Earthquakes `. ... 
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Theories of Scaling and Self- 
similarity 


George К. Zipf (1902-1950) 


C. 
and 

the | а 
to 


* 


my 


that Rank-order plot 
Zipf's law (1941) 
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Rank order of frequency 


review articles 


The language of genes 


David B. Searls 


Bioinformatics Division, Genetics Research, GlaxoSmithKline Pharmaceuticals, 709 Swedeland Road, PO Box 1539, King of Prussia, 

Pennsylvania 19406, USA (e-mail: david_b_searls@gsk.com) 

Linguistic metaphors have been woven into the fabric 140 
determination of the human genome sequence has brı 
imagination, with the natural extension of the notion ( 
life’. But do these analogies go deeper and, if so, can 
applied to molecular biology? In fact, many technique 
independently, may be seen to be grounded in linguist 
instrumental in extending our understanding of the Ізі 


— Protein kinase domain 
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100 E Sugar transporter domain 
WD domain 


__—— Helicase conserved domain 


he science of linguistics has fully as many facet 
ex ps and fields as biology, and like biology, what ma 
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be called its ‘modern era’ can be traced to th RNA-recognition motif 
1950s!. The decade that unveiled the structure c 
DNA also witnessed a revolution in linguistic 


led by Noam Chomsky, whose work radically diversifie: 


40 


Occurrence of words in Romeo and Juliet (+) 


20 


Occurrence of Pfam protein domains іп 5. cerevisiae genome &) 


Rank order of frequency 


Figure 3 Distributions of the number of occurrences of Pfam protein domains (blue 
squares) in the genome of the yeast Saccharomyces cerevisiae, and of words (red 
diamonds) in Shakespeare's Romeo and Juliet, in both cases sorted in rank order from 
left to right. The most frequently occurring domains and words are labelled. In both 
cases (and in many other genomes and texts) the curves are good fits to a power-law 
distribution known as Zipf's law, which relates the frequency to the inverse of the rank. 


Log-log plot 


A power law relationship between two variables 
yields 


a linear relationship between the logarithms of the 
quantities. 
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Zipf's law 


Zipf's law 
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log(rank) 
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* Esperanto 
* Latin 

* Ukrainian 
- Czech 

* Italian 


Spanish 
Slovene 
Finnish 
Hebrew 
Turkish 
Hungarian 
Galician 
Danish 
Belarusian 
Portuguese 


German 
Malay 
English 
Slovak 
Romanian 
Polish 
Uzbek 
French 


- Basque 

* Serbian 

* Dutch 
Catalan 

* Indonesian 
* Lithuanian 


Croatian 


City size distribution 


QUARTERLY JOURNAL OF ECONOMICS 
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FIGURE I 
Log Size versus Log Rank of the 135 largest U. S. Metropolitan Areas in 1991 
Source: Statistical Abstract of the United States [1993]. 
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Long Tail 


Why the Future of Business 
Is Selling Less of More 


Products 
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Rhapsody music downloads (for number of titles equivalent to Wal-Mart inventory) But it keeps going... (rank 25,000 to 100,000) 


180,000 
160,000 
140,000 
120,000 
100,000 
80,000 
60,000 
40,000 
20,000 


ы | | | | | 
0 5,000 10,000 15,000 20,00 25,000 


Downloads 
Downloads 


Track rank (thousands) 


Track rank 


and going (rank 100,000 to 800,000) 


Downloads 


Track rank (thousands) 


Force 1: Democratize the tools of production Force 2: Democratize the tools of distribution 


Result: More access to niches, 
which fattens the Tail 


Result: More stuff, which 
lengthens the Tail 


Force 3: Connect supply and demand 


Result: Drives business 
from hits to niches 


The evolution of the 80/20 rule 


Bricks-and-mortar retailer Long Tail retailer 


Products Revenues Profits Products Revenues Profits 


20% 


online-only 
inventory 


Vilfredo Pareto (1848- 
1923) 


Pareto principle 


Distribution of world GDP, 
198918! 


Quintile of population | Income 


Richest 20% 82.70% 
Second 20% 11.75% 
Third 20% 2.30% 
Fourth 20% 1.85% 
Poorest 20% 1.40% 


80/20 rule 


® Italian land ownership 
® British tax records 


$ Peapods in Pareto's 
garden 


Pareto distribution 


& Tail distribution 
® Power law tail 


® Long tail 


® Fat tail 
S a= 1+ 8 


$ epsilon is a symbol of a small number 


Тһе 80/20 rule 
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Тһе 80/20 rule 


Inequality and Gini coefficient 
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Gini = A/(A+B) | 


Cumulative % of X 


Normal distribution 


34.1%| 34.196 


Complementary cumulative 
distribution or Tail Distribution 


Ihe Power Law 
Distribution 


Power law distributions are ubiquitous, occurring in diverse 
phenomena, including city sizes, incomes, word frequencies, and 
earthquake magnitudes. 


Power laws easy to spot on log-log (doubly logarithmic) plots: 


Power Law PDF - Linear Power Law PDF - Log-Log 


Scale invariance 
and self-similarity 


Scale invariance 


® Pareto: distribution of wealth exceeding a minimum Xm 


® Change the SCALE and take those, who exceed x; 


& Pareto remains the same on the new SCALE Χ/Χι 


& There is a lack of natural SCALE in the problem 


& Only Power Laws are invariant under change of scale 


Scaling (scale 
invariance) 


& Distribution of height of people is NOT scale invariant. 

& Humans do have a typical height scale (1-2 meters). 

& There are no 10 m tall humans. 

& Tallness is not power law distributed 

& Number of people you can relate to (Dunbar number) 150 
& Number of years you can live (no one lived for 150 years) 
& Number of followers 1,10,100,1000,...,10 million 

& Number of sexual partners 1,10,100,1000 ... 


& Attraction, hype, fanship ... 


The Matthew Effect 


For whoever has will be given more, and they will have an abundance. Whoever does not have, even what they have will be taken from them. 


— 25:29, New International Version. 


The rich get richer and the poor get poorer 


& Robert K. Merton (1968) 
® Herbert A. Simon (1955) Nobel M. Prize (1978) 
& Udny Yule (1925) 


& Yule-Simon process and distribution 


Preferential attachment 


ko — 1: initial number of coins in box 
m: number of coins distributed in each round 
ko/m = «в can be a small number 


Small ball 
made of putty 


of different colours 
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(again, 2 times) 


FOLD & SQUASH 
(2 times) 
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(again, 2 times) 
ROLL UP 
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(again, 2 times) е 


FOLD & SQUAS CUT IN HALF 


(again, 2 times) 
~ FOLD & SQUAS 
(again, 2 times) 


Smale Horseshue 


Chaos, Mixing апа Entropy 


® Ludwig Boltzmann 

® Rudolf Clausius 

& Order and Chaos 

® 2nd law of thermodynamics 


& Entropy 


Тһе idea behind entropy 


During the last 150 
years ... 


& Statistical mechanics/physics was very successful 


® Laws of physics do not distinguish between past 
and future 


& Yet, we see that past and future are different 
& Things go from order to disorder 


& Coffee cools, sugar dissolves ... 


Тһе Reductionist Method 


& 5 July 1686 
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The first 300 years ... 


& Newton introduced the scientific method 
& Cut the system to elementary parts 

& These parts are featureless 

& Find out how the parts interact 

$ Newton's theory of gravitation 


® Elementary particles: electrons, photons, neutrinos, 
mesons, gluons, W and Z bosons 


® Very successful in the last 300 years 


Standard Model of 


FUNDAMENTAL PARTICLES AND INTERACTIONS 


The Standard Model summarizes the current knowledge in Particle Physics. It is the quantum theory that includes the theory of strong interactions (quantum chromodynamics or QCD) and the unified 
theory of weak and electromagnetic interactions (electroweak). Gravity is included on this chart because it is one of the fundamental interactions even though not part of the "Standard Model." 


matter constituents 


FERMIONS 


spin - 1/2, 3/2, 5/2, ... 
Leptons spin = 1/2 Quarks spin = 1/2 
5 Approx. ^ 
Flavor EEE Electric Flavor Mass Electric 
GeV/c^ charge сеу/с2 charge 


р. electron 
neutrino 


«1x10-8 


2/3 


€ electron | 0.000511 -1/3 Nucleus 


muon Size = 10-14 т 


U up 0.003 
d down 0.006 
nelle ine «0.0002 C charm 1.3 


2/3 


e 


|А. muon 0.106 S strange 0.1 -1/3 


Ln «0.02 t top 175 2/3 
neutrino 


T tau 1.7771 Б bottom 4.3 -1/3 


Spin is the intrinsic angular momentum of particles. Spin is given in units of h, which is the 
quantum unit of angular momentum, where fi = h/2x = 6.58х10-25 GeV s = 1.05х10-34 J s. 


Electric charges are given in units of the proton's charge. In SI units the electric charge of 
the proton is 1.60х10-19 coulombs. 


The energy unit of particle physics is the electronvolt (eV), the energy gained by one elec- 
tron in crossing a potential difference of one volt. Masses are given in GeV/c? (remember 
E = mc?), where 1 GeV = 10? eV = 1.60χ10-10 joule. The mass of the proton is 0.938 GeV/c? 
= 1.67х10-27 kg. 


Size < 10-19 m 


If the protons and neutrons in this picture were 10 cm across, 
then the quarks and electrons would be less than 0.1 mm in. 
size and the entire atom would be about 10 km across. 


force carriers 
spin = 0, 1,2, ... 


BOSONS 


Unified Electroweak spin - 1 Strong (color) spin - 1 


Electric 
charge 


Mass 
GeV/c? 


Electric 
charge 


Mass 
GeV/c? 


FOREN 
photon 
w- 


Name Name 


Electron 
Size < 10-18 m 


-1 Color Charge 

Each quark carries one of three types of 
“strong charge,” also called “color charge.” 
These charges have nothing to do with the 
colors of visible light. There are eight possible 
types of color charge for gluons. Just as electri- 
cally-charged particles interact by exchanging photons, in strong interactions color-charged par- 
ticles interact by exchanging gluons. Leptons, photons, and W and Z bosons have no strong 
interactions and hence no color charge. 


Neutron 
and 
Proton 

Size = 10-15 m Quarks Confined in Mesons and Baryons 

One cannot isolate quarks and gluons; they are confined in color-neutral particles called 
hadrons. This confinement (binding) results from multiple exchanges of gluons among the 
color-charged constituents. As color-charged particles (quarks and gluons) move apart, the ener- 
gy in the color-force field between them increases. This energy eventually is converted into addi- 
tional quark-antiquark pairs (see figure below). The quarks and antiquarks then combine into 
hadrons; these are the particles seen to emerge. Two types of hadrons have been observed in 
nature: mesons 44 and baryons 444. 


Residual Strong Interaction 

The strong binding of color-neutral protons and neutrons to form nuclei is due to residual 
strong interactions between their color-charged constituents. It is similar to the residual elec- 
trical interaction that binds electrically neutral atoms to form molecules. It can also be 
viewed as the exchange of mesons between the hadrons. 


PROPERTIES OF THE INTERACTIONS 


Baryons qqq and Antibaryons qqq 3 Mesons qq 
Baryons are fermionic hadrons. Interaction . . ! Mesons аге bosonic hadrons. 
There are about 120 types of baryons. Property Gravitationa There are about 140 types of mesons. 
Fundamental 
k Electri M: 6 5 i б 
Symbol Name content charge Gevie2 Spin Acts on: Mass - Energy Electric Charge | Color Charge | “irverscton Note Symbol Name content charge сеш SP 


Particles experiencing: 
proton 


1 Particles mediating: 
anti- 


proton Strength relative to electromag | 10718 m 


for two u quarks at: 


10-41 
10-41 
10-36 


3510 т 


for two protons in nucleus 


Matter and Antimatter 


For every particle type there is a corresponding antiparticle type, denot- 
ed by a bar over the particle symbol (unless + or - charge is shown). 
Particle and antiparticle have identical mass and spin but opposite 
charges. Some electrically neutral bosons (е.9., 20, y, and ης = cc, but not 
КО = ds) are their own antiparticles. 


Figures 

These diagrams are an artist's conception of physical processes. They are 
not exact and have no meaningful scale. Green shaded areas represent 
the cloud of gluons or the gluon field, and red lines the quark paths. 


A neutron decays to a proton, an electron, 
and an antineutrino via a virtual (mediating) 
W boson. This is neutron β decay. 


Graviton 
(not yet observed) 


Ап electron and positron *»— 


(antielectron) colliding at high energy can 
annihilate to produce B? and B® mesons 
via a virtual Z boson or a virtual photon. 


Quarks, Leptons Electrically charged Quarks, Gluons 
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60 
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Not applicable 
to quarks 
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The Particle Adventure 


A Visit the award-winning web feature The Particle Adventure at 
z0 http://ParticleAdventure.org 
This chart has been made possible by the generous support of: 
U.S. Department of Energy 
U.S. National Science Foundation 
Lawrence Berkeley National Laboratory 
Stanford Linear Accelerator Center 
American Physical Society, Division of Particles and Fields 


BURLE INDUSTRIES, INC. 


©2000 Contemporary Physics Education Project. CPEP is a non-profit organiza- 
tion of teachers, physicists, and educators. Send mail to: CPEP, MS 50-308, Lawrence 
Berkeley National Laboratory, Berkeley, CA, 94720. For information on charts, text 
materials, hands-on classroom activities, and workshops, see: 


http://CPEPweb.org 
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Two protons colliding at high energy can 
produce various hadrons plus very high mass 
particles such as Z bosons. Events such as this 


one are rare but can yield vital clues to the 
structure of matter. 


Emergence of Order from 
Chaos 


electrons and protons 


Why electrons апа protons 
can live forever? 


carbon synthesis 


why carbon exists іп the Universe? 


Emergence: Order from 
Chaos! 


&Formation of atoms from nucleons and electrons. 
&More and more complicated atoms... 
&Formation of molecules from atoms. 

Моге and more complicated molecules ... 
&Formation of condensed matter from molecules. 


ФМоге and more complicated forms of condensed 
matter ... 


®... 


&Formation of cells, tissues ... 


The science of complexity 


® Ilya Prigogine 


® Nobel Prize tn 
Chemistry (1977) 


® Entropy can decrease 
® Makes life possible 
® Order Out of Chaos 


In the last 40 years 


® | can increase my order if | can export my mess to 
you ... 


& Entropy inside can decrease if the system is open 
and can pump it out into the environment 


& Understanding open dissipative structures 


& Understanding the complexity of nonlinear 
dynamics, bifurcations, chaos. 


Evolution 


& Ok! We understand that it is physically possible to 
increase the complexity of systems 


& Ok! We understand how this happens technically 
& But, why the hell is this happening? 
& Any law of evolution? 


& Darwinian selection is a good start, but not enough 


The Santa Fe Institute 


Our researchers endeavor to understand and unify the underlying, shared patterns 
in complex physical, biological, social, cultural, technological, and even possible 
astrobiological worlds. Our global research network of scholars spans borders, 
departments, and disciplines, unifying curious minds steeped in rigorous logical, 
mathematical, and computational reasoning. As we reveal the unseen mechanisms 
and processes that shape these evolving worlds, we seek to use this understanding 


to promote the well-being of humankind and of life on earth. 


& Los Alamos National Lab 


& Civil utilization of research ideas of Cold 
War 


Phase transitions 


Ferromagnetism and the 
Ising model 
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Cluster size distribution 


10! 
Clustor sizo 


Mean field 
approximation 


Mean field diagram 


Bistability апа 
bifurcation 
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Epidemic processes 


Probability of propagation 


| 2 | | | 
9009 => өө | 0 τ 0.4 06 08 


Probability of infection 


SIS model 


Infective 
І 


Months with measles (%) 


Susceptible 
5 


Population size 
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